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Abstract 

An effective theory for the hole doped spin-1 antiferromagnetic chain is pro- 
posed in this paper. The two branches of low energy quasipaticle excitation 
is obtained by the diagrammic technique. In the large t limit (in which t is 
the hole hopping term), the lower band is essentially the bound state of one 
hole and one magnon and the other band is the sigle hole state. We find a 
critical value of t, t c = 0.21 Ah (in which Ah is the Haldane gap). For t > t c , 
with the decrement of t, the mixing of these two bands become stronger and 
stronger, and at the same time the effective band mass becomes larger and 
larger. When t < t c the minimum of the lower band moves away from zero the 
another point between zero and tt/2. The spin strcture factor is also calcu- 
lated in this paper, and we find that for large t limit the main contribution is 
from the inter-band transition which induce a resonant peak in the Haldane 
gap. While for small t limit the main contribution is from the intra-band 
transition which only cause a diffusion like broad bump in the Haldan gap. 
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I. INTRODUCTION 



Since Haldane's seminal paper on the quantum Heisenberg spin chains Q, it has been 
investigated by many theoretical and experimental studies. For spin-1 antiferromagnetic 
chain, it is now well established both experimentally and theoretically ||||] that there 
is an energy gap between the triplet excitation and the singlet ground state. The gap size 
is recently determined to be 0.41J by using density matrix renormalization group(DMRG) 
H and exact diagonalization ||. It is known that the low energy behavior for the quantum 
antiferromagnitc spin chains can be described by the nonlinear sigma (iVLcr)models and 
additional toplogical terms are further needed for half integer spin chains. Particular for the 
spin-1 chain, a phenomenological model named as free boson model has been proposed by 
Ian Affleck et al [p|-[T0|1 based on the large N expansion for the NLa model. It has been shown 
that the free boson model captures the basic physics of the spin-1 chain in low temperature 
and can be used to calculate many physical properties correctly. 

There are mainly two families of compounds which exhibit the essential physics of the 
antiferromagnetic quantum spin-1 chain. One is the early discovered compound CsNiCl^ 
and the other is Y2BaNiO^ . Replacing Y by Ca for Y2BaNiO§ |TTl , |T2ll one can dope 
the spin-1 chains by holes. The electronic transport properties, polarized x-ray absorption 
and neutron scattering of Y 2 - x Ca x BaNiO^ have been measured by J. F. DiTusa et al JT3 



The result of electronic resistivity shows that the hole doping greatly reduces the resistivity 
which implies a considerably large mobility for the holes. The neutron scattering experiment 
confirms the existence of new states in the Haldane gap fl3fl . Theoretically the effect of 
doping in Haldane gap system is a very interesting problem and have been studied by 
several groups using both numerical |14| and analytical methods |15|,|16 . 



For case of static hole doping, by using the DMRG method |14j Sorensen and Ian Affleck 
studied the in-gap state caused by two kinds of doping. One is bond doping and the other 
is site doping. In bond doping case, the effect of doping is to perturbe an antiferromagnetic 
bond with J' ^ J . It gives rise a localized state centered at the perturbed bond with a 
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descrete energy level appeared inside the Haldane gap for sufficiently strong or weak J'. 
In site doping case, a static hole is considered to be located on the O ion between two Ni 
ions. Therefore the super exchange process between these two spins is destroyed or partly 
destroyed by the hole. The bound states in the Haldane gap are found only when the 
coupling between the hole and the nearest spin is weaker than a critical value. Similar 



results were obtained by M.Kaburagi and et al using variational calculation |16 



Also one can consider a mobile hole doped in spin-1 chain. S.C. Zhang and D. P. Arovas 



17| considered a spin-0 hole hopping in a background of spin-1 chain by using an effective 
model which is quite similar to the t-J model and found some exact single and multi- hole 
spin singlet wave functions. The problem of spin 1/2 hole moving in spin-1 chain has been 



considered by K. Penc and H. Shiba |18| in the limit of small hopping amplitude. In their 



approach, the holes are hopping between the O sites and destroy the corresponding super 
exchange processes completely when the O sites are occupied by the holes. They found one 
spin-3/2 band and two spin 1/2 bands either in the VBS model or in the Heisenberg model 
due to the interaction between hole and its nearest neighboring spins. But for the realistic 
Heisenberg model, their variational calculation can only treat the finite lattice up to 15 sites. 
Recently E. Dagotto proposed a t-J like model to study the hole doped spin-1 chain. Using 
numerical techniques the dispersion of the effective hole bands and the spin structure 
factor were obtained for the finite lattice up to 12 sites. 

In the present paper, we propose an effective continuum theory to study the hole doping 
in spin-1 chain. The spin-1 chain is modeled by the free Boson model proposed by Ian 
Affleck et al |8|- [T0|1 which is essentially derived from the large N expansion of the NLa 



model. The effective interaction between the holes and the magnons are derived, based 
upon the following considerations. Recall the static hole limit, the effective interaction only 
contains the scattering process and the holes act as the scatters of the magnons. Then the 
effective Hamiltonian for the static hole doped in spin-1 chain can be even easily intepreted 
in the first quantizaction picture. The result shows that there are one or two (depends on 
the inetraction strength) impurity states with total spin equals 1/2 or 3/2 respectively in 



the Haldane gap which are actually the bound states of a hole and a spin-1 magnon. This 
result is in good agreement with the numerical results. We may then draw an intuition to 
the moving holes. When the hole moves,the effective interaction should contains two main 
terms, one describes the scattering process of the hole and the magnon which has the same 
origin with that of the static hole case, the other describes the emission and absorption 
process of magnons. 

We show in this paper that the first term results in a bound state of one hole and 
one magnon with the bound energy strongly depending on the total spin of the hole and 
magnon. Therefore, we will find that if we only consider the scattering process, the result 
is very similar to the static hole case. The only difference is that each impurity state in the 
static hole case will extend to a corresponding energy band if the hole can hop. The effect 
of the emission and asorption term will result a hybridization of single hole state and the 
bound state of one hole and one magnon with total spin 1/2. This term plays a crucial roll 
for small t, because the energy level of the bound state and that of the single hole state is 
very close in small t limit and is quite large in large t limit. Then we obtain three branches 
of excitation, the bound state of one hole and one magnon with total spin 3/2, the "one 
particle like" state with total spin 1/2 and the "two particle" like state with total spin 1/2. 
The "one particle like" state approaches to the single hole state in large t limit and the "two 
particle like" state approaches to the bound state of one hole and one magnon with total 
spin 1/2. The dispersion of these three branches of excitations are obtained in the whole 
range of the hopping amplitude t. The band minmum of the "two particle like" state is 
found to be located at 7r for t > t c and will be move toward 7r/2 when t < t c . The value of 
t c is found to be near 0.21A# in our calculation. In fact this can be understood as a result 
of the above mentioned hybridization effect in the small t case. These results are consistent 
with the results for VBS model obtained by K. Penc and H. Shiba. 

The spin structure factor is also obtained and the result is quite different for large t case 
and small t case. For large t case the hole contribution to the spin structure factor is mainly 
contributed from the interband transition(from "two particle like" state to "one particle 
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like" states) and the spectral weight whithin the Haldane gap is centered at the minimum 
energy cost between the two bands, which is consistent with the calculation of Dagotta et al. 
But for small t case the contribution is mainly from intra band transition(from "two particle 
like" state to "two particle like" state within the same branch of dispersion), so there exists 
a diffusion like broad bump in the Haldane gap. The difference of the spin structure factor 
in the above mentioned two limit is quite interesting and has never been mentioned in the 
previous works. 

The rest of the paper is organized as follows: The effective continuum Hamiltonian is 
derived in Sec. II, whereas dispersion of the bound states is claculated in Sec. Ill, and the 
spin structure factor is calculated in Sec. IV. Finally, we make a few concluding remarks in 
Sec. VI. 

II. The Derivation of The Effective Hamiltonain 

In this paper, we assume that the holes can hop whithin the oxygen sublattice, which is 
shown in Fig.l. 

Then we can begin with the following total Hamiltonian 

Hiatal = H ch + H h + Hj, + H Jl + Hj 2 

H c h is the rotationally invariant spin Hamiltonian for Antiferromagnetic spin-1 chain. 

H c h = J ^2 Si ■ S i+1 (la) 

i 

and H h is the hopping term of holes 

Hh = ~t E /„ + + i/ 2CT W,. + H.C. (lb) 

n 

Hji and H Jl represent the destroying of super exchange process and a Kondo like inter- 
action of spins and holes respectively, 

Hj> = - J' X] ■ 5 ? n+l/ T | + i/2, CT /n+l/2,<T (lc) 
<r,n 
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Hji — °n,<*,P ' (Sn + S n+ i) f^ n+1/2 f p, n+1 / 2 (Id) 

n,a,j3 

And Hj 2 describes the hole hop to another site with its spin flipped by interacting with the 
spins of the chain. 

Ej 2 = Ji &n,a,p ■ S n f+ n _y 2 fp >n+ i/ 2 + E.G. (le) 

n,a,/3 

It is known that E c h can be mapped to the non-linear a model by using the path 
intergral in the spin coherent state representation. The antiferromagnetic order parameter 
is represented by a three dimentional vector field H, obeying a constraint as |0| 2 = 1. 
The uniform magnetisation is represented by I, 

f=(J_U x ^ ( 2a ) 
V V dt v ' 

And the spin operator at site i can be written as, 

s(-iy$+t, (2b) 

Then the spin-1 Heisenberg Halitonian can be mapped into an effective continuum field 
theory with the Hamiltonian density as: 

**-? + §(g)' < 2c > 

with the constraint |0| 2 = 1. 

By taking the large N(the number of the components for field) limit, and further 
introduce a mass term as the Lagrange multiplier to relax the constraint on the field (j) into 
an averaged one, we follow the free Boson model or Ginsburg-Landau model proposed by 
I. Affleck et al |]8HlO|,in which an additional 4 term is added for keeping the stability. We 



then have 
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= i fl2 + I (I) + 1*** 2 + < 3 » 

where U.(x) is the canonic momentum conjugated to the field 4>{x) satisfying [<f>(x) a , U(x')/3} = 
^a,/3^{ x ~ x ')- The three parameters v , Ah and A are chosen phenomenologically to fit the 



experiment or the numerical result. In our present study we omit the 4 term as for a 
preliminary discussion. We can then diagonalize the Ti. c h by Fourier transformation, and 
obtain three branches of free magnons, which describe the triplet excitations upon the singlet 
ground state. 

^ = | E n(-*) • n(fc) + + ■ <fe 

k 

= Y. E k(<4^Kr + \) (4) 
k,r ~ 

In equation(4) dk, r a>t,r are the annihilation and creation operators of the magnons satisfying 

[a,k, r , at' A = Srr'Skk' 



with r — x,y, z, E(k) = \jv 2 k 2 + A 2 H and: 



= (afc + a-k) ( 5a ) 



fl(*) = <y§(3*-Sl fc ) (5b) 

Now we have expressed the Hamiltonian for the Heisenberg chain in terms of the three 
branches of gapful magnons, we have In the continuum limit, 

(Si + S i+1 ) w 2l( Xi ) + ^\x=xi (6) 

For the perfect spin-1 chain, the field l(x) is always one order smaller than the field 4>(x) 
and is of the same order as that of 9 tj^ which reflects the short range anti-ferromgnetic 
correlation in low temperature. But if doped with an spin- 1/2 holes, as we will show later, 
the hole will induce a localized mode of magnon which has an extension only of several 
lattice near the hole. So in the doped case near the site of the hole, the field l(x) could be 
of the same order with field <f(x), and could be one order larger than Therefore, we 

can omit in (6) and only keep the first term. Then the interaction term Hji, Hj 1 and 
Hj 2 can be written as, 



2 *r 



# * = 2 Ji E 



°a,l3f&i/2, a fi+l/2,l3 ~ 2^1 E ' ^a,/3/i+l/2,a/i+l/2,/J ( 7b ) 



(7c) 



We can now express the field /(xj), /(xj) 2 and 0(xj) in terms of the magnon creation and 
annihilation operators. We leave the detail derivation in appendix A. After discarding the 
multi-magnon processes, the Hamiltonian reads, 



H ch = E E k{at,n a k,n + 2- 

k,n 



(8) 



in which \i = —1, 0, 1 and 



a±{k) = [a x (k) + ± ia+(k) 



(9a) 



a o( k ) = a t( k ) 



(9b) 



Hj, = ~A E l'( k '^) J ' a tk'+q a ^'fa,k- q U,k (10) 

#Ji = E 2-f 1 (k',q)J 1 S^- Va,pa+ k , +q a U;k/ f+ k _ q fp, k (11) 

fiu,k,k' ,q,a,0 

as well as 

#J 2 = E 2 7i(fc', q)J2Cos(k)S llu ■ B a ,f)a+ iV+q a v ,k,f+ k _ q fp jk 

fiu,k,k' ,q,a,/3 

+ E 2J 2 sin(A;)^^ Qi ^(-g)/ a + (A; + g-vr)/ /3 (A;) + i/.a (12) 

in which D^ip =< jj,a\\j3 > is the transformation matrix between the (J 2 , J z ) representa- 
tion and that of (S z , a z ). 

In this paper, we choose the amplitude of the Haldane gap Ah as the unit of energy, 
so we have A H = 1. We choose the value of v by fitting the correlation length at zero 
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temperature, which is believed to be closed to 7 JT^|. So we have £ = = 7, then we 
have v = 7. The two vertex functions ji(k,q) and ^f'{k,q) are derived in appendix A. It is 
quite difficult for us to treat this k,q dependent vertex function analytically. So as a low 
energy effective theory, we simply replace the two k,q dependent vertex function by two k,q 
independent phenomenological parameters Ai and A' respectively, whose value can be fixed 
by fitting our approach to the existing numerical results for the static hole problem. 

In the static limit t = Ji = 0, the effective Hamiltonian is quite simple especially in the 
continuum limit, 



H = JdxJ2a;(x)(l- 



i d 2 



)a li (x) 



m b dx 

- 4 A' J' £ a+(0)a 7 (0) + £ 2AiJi£ 7At • s imp a+(0)a M (0) 



(13) 



in which 



dkan(k)e 



-ikx 



(14) 



To study the low energy excitation inside the Haldane gap, we need only consider the one 
magnon state. And this will lead to the following schrodinger equation in first quantized 
picture. 

i o 2 



2uih dx 2 



+ 2Ai JiSimp ■ S6(x) - 4X 2 J'5(x) + A H 



ip{x) = Eip(x) 



(15) 



The magnon wave function ip(x) has six components corresponding to six eigen states of sf mp 
and S z . We can easily solve the above schrodinger equation by first diagnolizing Hamiltonian 
in spin space. We have 

1 d 2 



+ 2J 1 X 1 e j 6{x) - A\ 2 J'5(x) + A H 



(16) 



2 m;, dx 2 

in which j is the total spin of the impurity and the magnon, and tj = j(j+l)/2— (S 2 +s 2 mp )/2, 



so 



ei/2 



-A 



H 



e 3/2 = 1/2 Ah 



(17) 
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The Eq.(14) describes a particle with mass m& bounded by a 5 potential which can be solved 
easily. The corresponding eigenvalues and eigen functions are: 

E 1/2 = A H - 2m b (A! Ji + 2\ 2 J'f i>y 2 {x) = -= exp~ N/L with L ' 



y/L 2m b {\ l J l + 2\ 2 J') 

(18a) 

E 3 / 2 = A - 2m 6 (-0.5Ai Ji + 2\ 2 J') 2 ikpix) = exp~ WL with L ' 



v 7 ! ^ m 6 (-AiJi + 4A 2 J') 

(18b) 

Eq.(18.b) is only meaningful for case of Ai Ji < 4A 2 J' If AiJi > 4A 2 J' the effective potential 
for the state with total spin 3/2 is repulsive which can be easily verified from eq.(16). 
Therefore there is no bound state with total spin 3/2 in this case. When AiJi < 4A 2 J' 
there always exist two bound state where one has two-fold degeneracy and the other has 
four-fold degeneracy. We can compare our results with the corresponding numerical results. 
Firstly we choose J' = J = 2.5 Ah which corresponds to the case that the super exchange 
interaction got destroyed entirely. If we have further J\ = 0, the problem becomes precisely 
equivelent to an open chain. The numerical studies show that the ground state of the open 
spin-1 chain is four-fold |2T||2"T| ] . In the simple approach engaged in this paper, it gives a 
nice description of these four-fold states. If we simply choose A 2 = 1.01 the energy cost 
to create one localized magnon near the chain edge is determined by (18. b), then the 
degenerated ground state will be zero magnon state as well as three- fold states with one 
localized magnon near the edge, which is altogether four-fold. The edge states decribed by 
(18. b) extend to 5 lattice which is also consistent with the numerical result. For finite J\, 
equation (18. a) leads to E\/ 2 ~ — 0.4AJi, E 3 / 2 m 0.2AJi in small J\ limit while the numerical 
study predicts that Ey 2 ~ — Ji, E$/ 2 ~ 0.5Ji. We then choose Ai = 2.5 to fit the numerical 
results. Fig. 2 shows the —E\/ 2 as the function of Ji in a full range of J\. Compared with 
the results obtained by Ian Affleck et al ||14|| , we find that our simple treatment fits quite 
well with the numerical results espcially in weak coupling regime. 
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III. The Dispersion of the Bound State of Magnon and Hole 



In this section we consider the bound state constituted by one magnon and one hole. 
In the present paper, we assume that the local super exchange is destroyed entirely, so we 
choose J 1 = J, Ai = 2.5, A' = 1.01, Ji = J 2 = 0.2A H , and m b = 1/49. We will consider this 
problem in two limiting cases, one is the large t limit and the other is small t limit. We can 
reorganize the total Hamiltonian as H = Hh + H c h + Hi + H 2 , in which Hh and H c h are the 
free Hamiltonian for holes and spins respectively as shown in (lb) and (4) and Hi and H 2 
are two kinds of interaction representing the scattering process and the magnon emission 
and absorption process respectively. 

Hi = ( 2A i J i + 2A i J2Cos{k)] ■ 

fiu,k,k' ,q,a,/3 

H 2 = £ g(k)D^ hp a+(-q)f+(k + q-n)Mk) + H.C. (20) 

in which g(k) = 2J 2 sin(fc) y^^f- 

a. The large t limit 

In latge t limit, the holes are mainly distributed near P = 0, so we can treat the bare 
hole dispersion in continuum approximation. We have E h (P) = —2t(cosP — 1) ~ m 
which to/ = Yf 

In order to study the dispersion relation of the bound state, we should consider the two- 
particle Green's function with total spin j(j=l/2 or 3/2) and its z-component m which is 
defined as the following. 

T jm (x,t) = -i ]T < jm\a'n' >< Ta^{x, t)f a ,(x, t)\a+(0, 0)/+(0, 0) >< afx\jm > (21) 

a'/i'afi 

The above two-particle Green's function can be obtained by using Fenyman diagram expan- 
sion. One may easily verify that the magnon emission and absorption term only act on the 
total spin-1/2 state and can't affect the state with total spin 3/2. So the spin-3/2 state only 
feels the scattering term as shown in Fig.3, which leads to, 
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So we have 



Tz(P,zv) = n (P,iu) + n (P,iu)T3(P,iu) (22) 



L3{P, IV) ~ 



1 - ir (P,iis) 



rn 1 • ////i'-p • 1 ) • /'/••(- -/ 

UPM= dqV^P-q) , {2 ™ b _ q)2 q2 (23) 



and 



to(p, -) = y ^ .. t , (P- g)a q2 A ; ( 24 ) 

in which V|(fc) = Ai [Ji + J 2 cos(A;)] — 4A 2 J' and P is the totale momentum for the hole- 
magnon system. 

In above equations ub and % are the Boson and Fermi distribution function respectively. 
But in our present study if we consider only the dilute limit of the hole doping, so we can 
assume that the temperature satisfy Tf <T<C A#, in which Tp is the Fermi temperature 
of the holes. Then we can replace the Fermi-Dirac distribution function n F (E) by its classic 
limit which is the Boltzman distribution function. In equation(26), we can further ignore 
the ub and zip in dilute doping and low temperature case. Then we have 

tiq(P, to ~\~ zO ) ~ dq 2 , 2 5 = / rfg ^ n 

J u + W-^-£- b -Delta H J-* q u, + zO+-^-£-A H 

(25) 

in which q' = q — j^-P and ^ = ^ + We may expand the functiom V| with respect 
to g' and only keep the leading order, because the main contribution to the integration 
comes from q' < ^2m b (A H + — a;), which is quite small in the a; regime in which we are 
interested. So we have 



Ren (P,u; + zO + )^--Vs[(l-^)P] i 



(26) 



^ Delta H - w + 



and 
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Reiio(P,u) = — - 



(27) 



for < A//. 

The bound energy for the bound state of one hole and one magnon with total spin-3/2 
is determined by the pole of T3(P,tu), which is 1 — Reno(P,uj) = 0. We then have 



E 3 / 2 (P) = 1 - 2^3/2 



+ 



2M 



(28) 



The situation of the bound state with total spin- 1/2 is not as simple as the case of 
total spin-3/2, because the interaction term Hj 2 will mix the two particle state and the one 
particle state, which is illustrated in Fig. 3. Therefore we should take into consideration the 
contribution from both the two terms. Following the diagrammatic rule shown in Fig. (4b), 
we then have 



r(p, w) = r (P, w) + iyp, w)G(p + n, i^r^p, %u)+ 

ri(P, iv)G{P + 7T, iu)T 2 (P, iv)G{P + 7T, ii/)ri(P, iv) + ... 



where 



r (P» = 



Ti(P,iv) = 



7T (P, iv) 



Vi 

2 



g(P,^P)n (P,iv) 



1 - Vi 



1-^)P 7T (P> 



and 



T 2 (P,tv)=g(P,^P) 2 7r (P,tv) + 



g(P,^Py4(P,iv] 


2 




a- 


"6 ^ p 

M ^ . 


1 - Vk 

2 






7T (P, 


iv) 



(29) 



(30) 



(31) 



(32) 



in which Vi(k) = 2X 1 [J\ + J 2 cos(k)] — AX' J'. To derive the above equation, the same 
approximation is done as (26). Then we have 

r (i» 



r(/» 



l-^(pap) G (p + 7r>a; )r (P,u;) 



(33) 



13 



in which 

G(P + ir,u) = (34) 

v ' ; lo - At + P 2 /2m f v ; 

is the Green's function near ir. Then we can obtain the energy of the bound state with total 
spin-1/2 by solving the eqaution: 

1 - g 2 (P, ^P)G{P + tt, w )r (P, a;) = (35) 



The result is 



i P 2 

Ei (P) = A H - -fiV? + — — (36) 
2 V ' 2 2 2M* 

with the renormalized effective mass 

1 1 3J|(l-|&) 2 &iVi(0) 

in which m/ = ^. We find that if t is smaller than a critical value t c , the effective mass will 
change sign which means the band minimum may move away from the P = state and will 
be located at a point between and n. From the above expression, the t c is determined by 
the condition = 0. And the result is t c = 0.21A#. 

For the case of t » t c , the magnon emission and absorption processes are no longer 
important and can be ignored, so the number of the magnons is conserved. Then the essential 
physics can be viewed more clearly as a two-body problem in the first quantization picture. 
We can write again a Schrodinger equation for the hole-magnon system. 
Id 2 Id 2 



2m b dx 2 2rrifdxf' 



+ Js imp ■ S5(x b - Xf) - 4X 2 J'S(x b - Xf) + 1 



i/j(x b ,x f ) = Ei/j(x b: x f ) 
(38) 



in which mj = ^ and J = 2Ai(Ji + J2). The above Halmiltonian describes a two body 
problem with an attractive interaction between them. We can divide the Hamiltonian into 
two parts, one describe the motion of center of mass and the other describe the relative 
motion. 
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H r = + Jsim P ■ SS(x) - A\ 2 J'5(x) + 1 (39b) 

with // = ™ f J^ b ,M = m b + rrif,x = Xf — x b and x = mbXb ~^ n f x f m Equation (23a) describes 
a free propogation of a composite particle constituted by a hole and a magnon with total 
mass M, and equation (23b) describes a particle in a attractive potential which is solved in 
(19a) and (19b). So the situation in large t limit is quite clear. One magnon and one hole 
form a bound state with the bound energy 1 — 2/i(— 0.5Ai J\ + 2X2J') 2 for S to tai = 3/2 state 
or a bound energy 1 — 2//(Ai Ji + 2A 2 J') 2 f° r Stotai = 1/2 state. These results are consistent 
with and further confirm the results obtained by using Feynman diagram techeniques in the 
case of t » t c . 

The center of mass moves like a free particle with total mass M. So in the case of large 
t limit, the bound state of magnon and hole with total spin 1/2 is energetically favorable 
against the free hole state. Therefore, the ground state of the system will be the bound pairs 
of the holes and magnons with total spin 1/2, in another words, every hole will induce a 
magnon binding with it. Further, these bound pairs can be viewed as a composite Fermions 
with spin 1/2 which can propogate freely in the spin-1 chain. 

b. The Small t limit 

In the small t limit, we should calculate the same diagrams as in large t limit. Also we 
can calculate the dressed two particle Green's function by considering the diagramm shown 
in Fig(4). The difference between the two limiting cases is that we can't use the long wave 
length expansion anymore for holes in the small t case. The it (P,uj + 10 + ) in small t limit 
can be written as: 

n (P,uj + iO + ) = I dq 5 (40) 

J u + tO+ + {2tcos(P-q)-2t)-£- b -A H V ; 
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Since in small t limit we always have l/m& ^> t, the main contribution for the integrand 
comes from a quite narrow q regime near q equals to zero. So we can expand cos(P — q) in 
q near P, that is 

1 



• /(I' ) Jdq , , . ,, , ,, ., „, 



uo + i0+ + 2t cos(P) + 2t sin(P) -q-t cos(P)g 2 - ^ - A H - 2t 

(41) 

then intergrate it using the method described before. Also we can obtain the values of ~K\ and 
7r 2 using the same method. The dispersion relation of the qusiparticle which is the mixed 
state of the single hole state and the bound state of one magnon and one hole is determined 
by searching for the low energy poles of the Green's function similar to what we have done 
in the large t limit. Two branches of excitation are obtained with energy u±(P) as shown 
in Fig. 5 for various value of t. Near the poles, the two particle propagator can be written as 

T(P,lu + iO + ) = Zi ^ P) - (42) 

for « = +,—. The result of Zi(P) is shown in Fig.6. From Fig.6 we can find that for P 
near zero the lower branch is almost two particle like and the higher branch is one particle 
like, but for P near 7r the lower branch is one particle like and the higher branch is two 
particle like. Between the above two limit, the one particle state and two particle state 
are strongly hybridized. So near and n, the mixing of the one particle state and the two 
particle state is quite small and the effect of the magnon emission and absorption term is 
only modifying the effective mass of the bounded states. But for P near tt/2, the one particle 
state and the two particle state are strongly hybridized, which may split the two energy level 
remarkably. When t is large this hybridization effect only increases the effective mass. But 
if t is sufficiently small the energy split caused by J 2 term which is of the order J 2 could 
be much large than the kinetic energy of order t. Then the band minimum of the lower 
excitation will no longer stay at P = but move towards n/2. We can find from Fig.5 the 
critical value of t is 0.21A#, which is consistent with the value that we obtained in the large 
t limit calculation. From equ.(12) we can see that the origin of the magnon emission and 
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absorption term is the coupling between the mobile holes and the antiferromagnetic field 
4>(x). Although the long range antiferromagnetic order is not stable in ID spin chain, the 
short range antiferromagnetic fluctation is still quite strong in low temperature. If the holes 
can feel such antiferromagnetic fluctation, the band minimum will move to n/2. It is rather 
interesting to note that this feature is quite similar to that of the single hole moving in the 
2D antiferromagnetic lattice. 

Also, we can calculate the renormalized single particle Green's function considering the 
diagrams showing in Fig(4c). And the same quasiparticle poles are obtained. But we provide 
here a transparent intuitive picture. 

Our results for small t limit are very similar with the results obtained by Shiba el al |18| 



for the VBS model, in which the band minimum is moved towards 7r/2 when the hopping 
term t is smaller than a critical value. 



IV. Calculation of the Spin Structure Factor 

The spin structure factor can also be calculated in our approach. Since in the doped 
spin-1 chain both the holes and the spins can contribute to the spin structure factor, so we 
have 

k k 

(43) 

So the spin structure factor can be divided into three parts, the magnon part S m , the hole 
paet Sh and the mixed part. We can prove that the contribution from the mixed term is zero, 
because the magnon emission and absorption vertex function satisfy g{k, 0) = —g(—k,Q). 
Using the fluctuation-disipative theorem, we have 

S m {-K,u) = 2>[ImD x {ti,u) - ImD x {p,-u)\ [1 +n b (u)] (44) 

in which Z?i(0, u) is the dressed magnon Green's function with S z = 1. 



17 



S h (n,uj) = Imxh(ir,u})n b (uj) 



(45) 



in which Xh(j, uj) is the spin susceptibility of holes. 

Firstly the dressed Green's function of magnons is calculated by considering the self 
energy diagram in Fig.7(a,b). Fig. 7(a) represents the contribution of scattering term and 
the A S (P, uj) is the effective scattering potential of the holes and magnons, which can be easily 
derived from the full two particle Green's function T S (P, uj) as A s (P,uj) — V s + V s T(P,uj)V s 
in which index s represent two different channels with total spin 1/2 and 3/2. The Fig. 7(b) 
represent the four kinds of contribution from the magnon emission and absorption term 
which correspond to single particle to single particle, single particle to two particle, two 
particle to single particle and two particle to two particle transition respectively. And the 
diagrams considered by us to calculate Xh(j, uj) are shown is Fig.7(c), which is quite similar 
with Fig. 7(b) except that there are no vertex function g(k). 

Since in this paper we are only interested in the low energy response within the Haldane 
gap, we can discard the high energy part during our calculation. When t is much larger 
than t c , the band minimum of the lowest excitation is at zero. Therefore the diagrams in 
Fig.7(b,c) describe the intra band transition from zero to ir which costs the energy 4t. In 
the large t limit this energy scale is much larger than the Haldane gap, so the intra band 
transition can be ignored in the large t limit. Therefore the only important diagram in large 
t limit is Fig. 7(a), which is the inter band transition from the two particle bound state to 
the one particle state. According to Fig. 7(a) the self energy of the magnon can be written 
as, 




a,m,s,n 



(46) 



in which [i and a are the spin index of magnon and hole respectively. 



A s (iz/ n , P) can be expressed by its spectral function as 




(47) 
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in the above expression the first term is contributed by the bound state with A s = [iV^ 
and Pciyj') and the second term which is nonzero above the Haldane gap is attributed to 
the scattering state of one magnon and one hole. Since we are only interested in the spin 
response whithin the Haldane gap, we can omit the continuous part and only keep the first 
term in the above equation. Then we have 

ImE(u,iJi,k)= [dPj2 ~KD% 

a,m,s 

^- A « +1 2" V - + & t -^ ] <48) 
The spin response function whithin the Haldane gap for t=2 (large t limit) is shown in 
fig.8(a), we can see clearly that there is a resonate peak in the Haldane gap. This low energy 
peak can be attributed to the transition from the bound state of one magnon and one hole 
with total spin 1/2 (which is the ground state under the parameters we chosen here) to the 
free hole state. 

The situation of small t limit is quite different with the large t limit. Now the band 
minimum of the lower state is moved toward ir/2. So unlike the large t limit, the holes dressed 
with magnons are distributed near ir/2. Therefore the intra band transition described by 
fig. (7b) and fig. (7c) become very important now, because the momentum transfer n costs 
very small energy for the states near tt/2. While the energy cost for the inter band transition 
is increased by the virtual magnon emission and absorption process. So in small t limit all of 
the diagrams in fig. 7 must be considered. We have calculated all these diagrams numerically, 
and the result is shown in fig. (8b) for t = 0.1 Ah- The contribution from all these figures 
are evaluated under the same approximation as the large t limit which keeps only the low 
energy peak in the spectral function. For example the contribution from the second diagram 
in fig. 7(b) can be written as 

Ei(iw,7r)= / > rfP^7 2 (P,iz/,^)Ai(P,iz/)Ai(P + vr,iz/ + ^) (49) 



AP-k)\ 1 2 P 2 , 
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in which 



l(P, ™y iuJ ) = I 7T-3 dk Z " ~77T ' " — Ti i — V ' ~ : \ (P^kW 

J 2vr/3 *vi ~ e{k) ito + iv 1 - e{k + vr) w _ - A H - 



dA; 1 1 

-{n F [e(A ; )]^ + e(A;) _ e(fc + 7r) ■ tu _ e{k) _ AH _(P^E 

1 1 

+n F [e(fc + 7r)] 



+ 



1 + n b (A H + 



e(k + tc) -iu- e(k) % v + i u - e (k + tt) - A H - 
{P-kf 



iu-A H -^-e(k) ^v + ^u J -e(k + n)-A H -^ J 

(50) 



whith e(k) = — 2t[cos(k) — 2]. The first two terms in the above equation can be ignored if 
we only consider the low temperature case. Also we use the same approximation that only 
the low energy peak of X(P, iu) is considered. Then after the summation of iu has been done, 
we have 

ImZtiu,*) = f dPj2S(co-cu m (P + n)+uj n (P))A n {P)A m (P + 7r) 

mn 

[n F {u n {P)) - n F {u rn {P + vr))] 7 2 (P, u n {P),u) (51) 

in which u n (P) is the dispersion relation of the two branches obtained in equation(42). The 
other diagrams in fig. 7(b) can be calculated similarly. Finally we calculated the numerical 
result for the spin response function in small t limit. The result is shown in Fig. 8(b). We 
can see clearly from fig. 8 that a diffusion like peak is present at low energy near zero resulted 
from the intra band transition in small t limit. The energy of the inter band transition is 
near 1.22 Ah which is out of the range in which we are interested. 

Our results in large t limit are consistent with the results of Dagoto and et al which have 
two peaks represent the bound state to single hole state and the intrinsic magnon excitation 



respectively. In the paper of Dagoto and et al []TS] , the dynamically spin structure factor in 
small t limit has not been calculated. So our calculation is the first work which indicate the 
difference in the dynamically spin structure factor for large t limit and small t limit. We 
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find that when t is smaller than a critical value, the main contribution to the dynamically 
spin structure factor will change from the inter band transition to intra band transition. 
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Appendix 

In the appendix, we derive the representation of field l(q) and l 2 (q). 
First we have 



-> -> -> i IE 

Kq) = + s) x n(-fc) = E 2\ ^^^ k+q + ^ k -^ x ~ ^ 



= S W Y^(- d k+ q xat + at k _ q x a_ fe + 4 +g x a_ fc - ^ fc _ 9 x 3+) (52) 

It is quite clear that the first two terms describe the scattering process and the last two terms 
describe the multi-magnon processes which is unimportant when we are only interested in 
the low energy physics. So the terms describing the multi-magnon processes can be omitted 
in our present paper. And we have 

Kq) ~ J2li( k ^) s t x a k+q (53) 



k 



in which 71 (k, q) — i 



k+q 



— ^2 7i(&'> k +Q)li(k", — &)•[— 2a+ jfc ,a yife »_ fc 5 fe / +fe+ q ife »— 2a+^ 

U UI h.II 
ru 1 re j re 



Jra t,k' a t,k" a y,k'+k+ q a y ,k" -k + 0'y^'0-y !k iiClz, k '+ k +qO-z, k "- k + ^t,k' a t,k" a x,k'+k+q^x,k"-k 



+0'y,k' a y,k" a x,k'+k+qO'X,k"~k + & z,k' a z,k" a y,k' +k+qO-y,k" -k + CL Xik ' a x,k" a z,k'+k+qO-z,k"-k 



' a t,k' a t,k" a x,k'+k+qa y ,k"-k — at, k ' a y, k " a y,k'+k+qaz,k"-k — at, k ' a t,k" a z,k'+k+ q ax,k"-k 
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~ a t,k' a y,k" a y,k'+k+ q ax,k"-k — a y,k' a t,k" a z,k'+k+ q a y ,k" ~k ~ a-t,k' a t,k" a x,k'+k+q^z,k"-k] 

After discarding the multi-magnon terms, we have 

l 2 (q) = Y,2l'(k',q)al kf a a , k , +q (54) 

k',ce 

with 7' = Ekli(k' ,k + q^k' + k + q, -k). 
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Figure Caption 



Fig.l The illustration of holes doped antiferro magnetic chain. 

Fig. 2 The camparison of our result and the numerical result for the lowest excitation in 
static hole limit. The line represents our result and the squares represent the corresponding 
numerical result. 

Fig. 3 The Feyman diagrams used in this paper, (a) The full line represents the hole's 
Green Function. (b)The dashed line represents the magnon's Green function, (c) The 
scattering vortex contributed by H ± . (d) The magnon emission and absorption vortex con- 
tributed by Hi. 

Fig. 4 (a) The dressed two-particle Green function with total spin 3/2. (b) The dressed 
two-particle Green function with total spin 1/2. (c) the dressed hole's Green function. 

Fig. 5 From down to up, the up three curves represent the uj + (P) for t=0. 22, 0.20,0. 18 
respectively. And from up to down, the lower three curves represent the u>-(P) for 
t=0.22,0.20,0.18 respectively. 

Fig. 6 The full and dashed lines represent the spectral weight z+(P) and z-(P) respec- 
tively. 

Fig. 7 (a) The magnon self energy caused by the scattering term, (b) The magnon self 
energy caused by the magnon emission and absorption term, (c) The hole's contribution to 
the spin susceptibility. 

Fig. 8 (a) The spin structure factor in large t limit (t=2). (b)The spin structure factor 
in samll t limit (t=0.1). 
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